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Abstract 

In this paper, we consider spatially homogenous and anisotropic 
Bianchi type I universe in the context of F(T) gravity. We construct 
some corresponding models using conservation equation and equation 
of state parameter representing different phases of the universe. In 
particular, we take matter dominated era, radiation dominated era, 
present dark energy phase and their combinations. It is found that 
one of the models has a constant solution which may correspond to the 
cosmological constant. We also derive equation of state parameter by 
using two well-known F(T) models and discuss cosmic acceleration. 
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1 Introduction 



Modified theories of gravity have recently gained a lot of interest due to its 
possible explanation about dark energy (DE). Modern cosmology is in a state 
of crises in a sense that it started with dark matter and erupted with an in- 
dication that most of the universe is made up of DE. The discovery of the ac- 
celerated expansion of the universe [l|-[3] indicates that the universe is nearly 
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spatially flat and consists of about 74% DE causing cosmic acceleration. This 
unknown energy (having negative pressure) is physically equivalent to vac- 
uum energy and is almost equally distributed in the universe. It has been 
used as an essential factor in a recent attempt to formulate a cyclic model 
of the universe. The universe (in its developing process) passes transiently 
through the stiff fluid era (u — 1), the radiation dominated era (u = 1/3), 
matter dominated era (u = 0), transition era (u = —1/3) and then tends to 
the DE dominated era (u = —1). In General Relativity (GR), the simplest 
and most appealing candidate of DE is the cosmological constant. However, 
it suffers from two serious theoretical problems, the cosmological constant 
problem and coincidence problem [I]. 

In alternative theories, f{R) theory has many applications in cosmology 
and gravity. This theory can directly be achieved by replacing the Ricci scalar 
R by f(R) in the Einstein-Hilbert action in GR. The study of the physics 
of f(R) models is, however, hampered by the complexity of the fourth order 
field equations in the framework of metric formalism [5]-[Tl]. Following the 
same scheme of modification in action, F(T) gravity can be obtained by 
replacing torsion scalar T by its general function F(T) in the Lagrangian 
of teleparallel gravity [15]-[21j. It helps to explain the accelerated universe 
without introducing any DE component |22j-[25j. 

The F(T) gravity models use the Weitzenbok connection which has no 
curvature but only torsion. Here torsion is responsible for the accelerated 
expansion of universe and is formed using four parallel vector fields, called 
vierbiens [26], which are linearly independent. In this framework, the torsion 
tensor is formed from the products of the first derivatives of tetrad. An 
important advantage of this theory is that its field equations are of second 
order and hence easy to tackle as compared to f{R) theory. 

Myrzakulov [27J discussed different F(T) models including scalar fields 
and gave analytical solutions for the scale factors and scalar fields. The same 
author studied the relationship between F(T) gravity and /c-essence [28J and 
also presented some new models of purely kinetic fc-essence. Karami and 
Abdolmaleki [2H] obtained equation of state (EoS) parameter of polytropic, 
standard, generalized and modified Chaplygin gas in this modified gravity 
scenario. Wu and Yu [30] discussed the two new F(T) models and showed 
how the crossing of phantom divide line takes place. They also discussed 
the observational constraints corresponding to these models. Dent et al. [3T] 
discussed F(T) cosmology both at background and perturbed level. They de- 
rived expressions for growth factor, stability of this theory and vector-tensor 
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perturbations. Li et al. [32] explored this modified gravity and local Lorentz 
invariance and remarked that this theory is not local Lorentz invariant. Yang 
[33] introduced some new F(T) models and described their physical impli- 
cations and cosmological behavior. Bengochea [M] investigated the conse- 
quences of data sets in this modified gravity. All the above mentioned work 
have been carried out for the FRW metric. 

In this paper, we would reconstruct the F(T) gravity models using Bianchi 
type / spacetime which is the generalization of FRW metric [35] . This the- 
ory becomes equivalent to GR if F(T) is replaced by a constant [3~T]-|32j. 
The paper is organized as follows: In next section, we present some basics 
of teleparallel gravity and the corresponding field equations for Bianchi / 
are given in section 3. A detailed construction of F(T) gravity models is 
given using two approaches in section 4. Section 5 is devoted to study the 
EoS parameter for two particular models and also a discussion on cosmic 
acceleration is provided. In the last section, we summarize and conclude the 
results. 



2 Preliminaries 

In this section, we introduce briefly the teleparallel theory of gravity and its 
generalization to F(T) theory. The Lagrangian density for teleparallel and 
F(T) gravity are, respectively, given as follows 



L ^> = ih F(T) - (2) 

where T is the torsion scalar, F(T) is a general differentiable function of 
torsion, G is the gravitational constant and e = \J—g. Mathematically, the 
torsion scalar is defined as 

T = S p ^ U T P p U , (3) 
where S p pv and torsion tensor T p ^ are given as follows 

S P ^ = ^{K pu p + 5 p T dv e — 5 u p T ep e ), (4) 
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Here h l are the components of the non-trivial tetrad field hi in the coordinate 
basis. It is an arbitrary choice to choose the tetrad field related to the metric 
tensor g pv by the following relation 

g»v = Vijh^K, (6) 

where rjij = diag(l, — 1, — 1, — 1) is the Minkowski metric for the tangent 
space. For a given metric there exist infinite different tetrad fields h 1 ^ which 
satisfy the following properties: 

Kh» = 6}, h^ = 5;. (7) 

The procedure for evaluating the tetrad field has been given in many papers 
[T5]-[2I]. Notice that the Latin alphabets (i, j, ... = 0, 1, 2, 3) will be used to 
denote the tangent space indices and the Greek alphabets (/z, v, ... = 0, 1, 2, 3) 
to denote the spacetime indices. The contorsion tensor K^ u p is defined as 

R ,u p = _ 1 (TM , p _ T „v p _ ^ ^ (g) 

which is equal to the difference between Weitzenbock and Levi-Civita con- 
nections. The variation of Eq.(|2]) with respect to the vierbein field leads to 
the following field equations 

[e-'d^eS, n - h^T" pX S p "»]F T + S t ^d p {T)F TT + \h\F = ~K 2 h?T v p . (9) 

Here F T =%, F TT = k 2 = 8irG, S 4 ^ = h?S p ^ with antisymmetric 
property and T pv is the energy-momentum tensor given as 

Tp = diag(p m , -p m , -p m , -p m ), (10) 

where p m is the density while p m is the pressure of matter inside the universe. 

3 The Field Equations 

The line element for a flat, homogenous and anisotropic Bianchi type / uni- 
verse is 

ds 2 = dt 2 - A 2 (t)dx 2 - B 2 {t)dy 2 - C 2 (t)dz 2 , (11) 
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where the scale factors A, B and C are functions of cosmic time t only. Using 
Eqs.flS]) an d (fTTT) . we obtain the tetrad components as follows 



h % = diag(l, A, B, C), /if = diag(l, A~ , B~ L , C 



(12) 



which obviously satisfies Eq.flT]). Substituting Eqs.(j3J) and (JSJ) in (j3J) and 
using pip , it follows after some manipulation 



T = -2 



/ AB BC CA 
I AB + B~C + CA 



(13) 



The field equations (Q for i — = v and i — 1 = v turn out to be 



, AB BC CA , 



(14) 



AB BC CA B C\ IB C 

U AB +2 BC + CA + B + c)^- i \B + C 



X 



A A 2 
A ~ A? 



B C 
B + C 



C C 2 \ I A B 
1 C ~ C 2 I \A + B 



B _B 2 
B~B~ 2 

F TT -F = 2k 2 Pt , 



C A 
C + A 



(15) 



The conservation equation takes the form 



. A B C . 

Pm + ( -J + "g + -q ) {Pm + Pm) = 0. 



(16) 



The average scale factor R, the mean Hubble parameter H and the 
anisotropy parameter A of the expansion respectively become 



R 
H 



(ABC) 



1/3 



If A B C\ _R 
3\A + B + C ) ~R' 



(17) 
(18) 



i=l 



(19) 



5 



where are the directional parameters in the direction x, y and z respec- 
tively given as 

Hl= A, ft=§, H 3 = ^. (20) 

It is mentioned here that the isotropic expansion of the universe is obtained 
for A = which further depends upon the values of unknown scale factors 
and parameters involved in the corresponding models [37j-[39j. Equation 
( TT3|) can be written as 



implying that 
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If we take F(T) = T, then Eqs.flHD and (USD will reduce to 



1 

Pm + PT= 2^ 



Pm+PT = 77Z2 



1 

2T 2 



, AB BC CA , 

-4 1 1 + T 

AB BC CA 



f AB BC CA B C 

2 h 2 1 1 1 

\AB BC CA B C 



T 



where pr, Pt are the torsion contributions given by 
1 



Pt 



Pt 



, AB BC CA 
-4 I — + — + 



+ 




(21) 
(22) 

(23) 
(24) 



(25) 



(26) 



The relationship between energy density p and pressure of matter p is 
described by EoS, p = up, where u is the EoS parameter. For normal, 
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relativistic and non-relativistic matters, EoS parameter has different corre- 
sponding values. Using Eqs.( IT3|) and (1T5]) . the EoS parameter is obtained as 
follows 

(AY - 2E)F T - 4ZF TT 



-1 + 



-AUF T + F 



where 



E 
Y 
Z 



AB BC CA 
+ 2— + 



AB 
BC 



BC CA 



B 
B 



C 
C' 



B C 
BC + B + C' 



(27) 

(28) 
(29) 



x 



U 




(30) 



(31) 



It is mentioned here that the homogenous part of Eq. (jl4p yields the following 
solution 

F(T) = (32) 



T 

where Cq is an integration constant. Using this equation in Eq. (fl5|) . we obtain 

1 ( 3MT 3H + J + L 1 



Vr, 



c 



2k 2 \ 2T 2 



T 



T 



(33) 



where L = — ^ and M = § + We would like to mention here that p Ti 
vanishes for the FRW spacetime 



4 Construction of Some F(T) Models 

Here we construct some F(T) models with different cases of perfect fluid by 
using two approaches. In the first approach, we use the continuity equation 
ffT6j) while in the second approach, EoS parameter fT27]) will be used. As the 
constituents of the universe are non-relativistic matter, radiations and DE, 
we consider the corresponding values of u in the following subsections. 
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4.1 Using Continuity Equation 

In this approach, we use the following relation jlQ] for the Bianchi type I 
universe 

9{a + B + c) - Ho + 3ABC> (34) 

where H is the Hubble constant having primary implication in cosmology 
and p is an integration constant. The value of H corresponds to the rate 
at which the universe is expanding today. This equation implies that 



(ABC)-i = ^-(H*-H*). (35) 



3 

« 2 Po 

Using EoS in Eq.f lTB]) . it follows that 

+ 3/7(1 + u) = 0. (36) 

Pm 

The components of the universe are described by the terms dark matter and 
DE. We consider different cases of fluids and their combination to construct 
corresponding F(T) models. For example, for relativistic matter, u = 1/3, 
for non-relativistic matter, it is zero and for DE era, it is equal to —1 [4"T] . 

Case 1 (u = 0): 

This is the case of non-relativistic matter, like cold dark matter (CDM) and 
baryons. It is well approximated as pressureless dust and called the matter 
dominated era. Inserting uj = in Eq.(|36l) and using ( |35l) . we have 



p m = p c (ABC)- 1 = 4r^(# 2 - #o 2 )> (37) 

where p c is an integration constant. In terms of torsion scalar, above equation 
becomes 

p m = 7 £-(J-9H*-T). (38) 



Substitution of this value of p m in (1T41) implies that 

2 Pc , T n Tj2 



2TFt + F = -—(J — 9Hq — T) (39) 
3po 
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which has the solution 



F(r) ^(^/^- iM »-H' (4o) 

This will have a unique solution if the value of J is known which corresponds 
to the unknown scale factors. Thus for matter dominated era, we obtain a 
model in the form of torsion scalar and Hubble constant. 



Case 2 (u = 1/3): 

Here we consider the relativistic matter, like photons and massless neutrinos, 
with EoS parameter u = 1/3. This case represents the radiation dominated 
era of the universe. Substituting cj = 1/3 in Eq.( l36|) and using ( 12T|) and (|35|) . 
we obtain 

where p r is another integration constant. Inserting this p m in ( IT4T) . we get 



implying that 



2TF T + F = ~^-j( J - 9H 2 - T) 4 (42) 



8lK<>p' VTj yp/" TVS 



This also depends upon the value of J as well as torsion scalar and Hubble 
constant. 



Case 3 (to = -1): 

This case represents the present DE constituting 74% of the universal density. 
Dark energy is assumed to have a large negative pressure in order to explain 
the observed acceleration of the universe. It is also termed as energy density 
of vacuum or cosmological constant A. Replacing u = — 1 in Eq. (!36|) . we 
have 

Pm = Pd, (44) 
where pa is an integration constant. Consequently, Eq.( TT4"]) takes the form 

2TF T + F = 2k 2 p d (45) 
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with 

F(T) = 2k 2 p d . (46) 

This turns out to be a constant model which is consistent with cosmological 
constant. 

Case 4 (Combination of to = and to = 1/3): 

Let us now consider the case when the energy density is a combination of 
two different fluids, the dust fluid and radiations. Adding Eqs. (l38|) and (j41j) . 
after simplification, it follows that 

*• = sk (J ~ ml ~ T)k + ^ky ~ 9H <> ~ T)% (47) 



Substituting this p m in dHJ), we get 

2TF T + F = A (J _ 9H * _ T)[p c + z^(J ~ 9H 2 - T) 3 ] (48) 
and hence 



F(T) = £±(J- [ jLdT - 18H 2 - -T 



3p \VT J VT 

+ Pi / Y— -HS-T^Y dT (49) 

+ SIk^VtJ \TV* TVs J ^ 

Case 5 (Combination of to = and w = — 1): 

The combination of EoS parameters for matter dominated era and DE yields 

Pm = tt4-( J - 9H 2 -T) + Pd . (50) 
3k 2 Pq 

Inserting in Eq. ffT4"]) . it follows that 

2TF T + F = p^(J - 9H 2 — T) + 2k 2 p d (51) 



yielding 



F{T) = t{7rS^ iT - l&Hl -\ T]+2 ^ (52) 
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Case 6 (Combination of uj = — 1 and oj = 1/3): 
This case gives the following form of the energy density 

p m = p d + ^- A {J-$Hl-T)\ (53) 
Substituting this value in Eq.flHJ), we get 

2TF T + F = 2k 2 p d + ~^-,{J - 9H 2 - T) 4 (54) 

which gives 

F(T) = 2k 2 p d + [ T l A 4 dT. (55) 

It is mentioned here that the cases 4-6 provide F(T) models for combi- 
nation of different matters. Normally, the dark matter and DE developed 
independently However, there are attempts |42] to include an interaction 
amongst them so that one can get some insights and see the combined effect 
of different fluids. Dark matter plays a central role in galaxy evolution and 
has measurable effects on the anisotropics observed in the cosmic microwave 
background. Although, matter made a larger fraction of total energy of the 
universe but its contribution would fall in the far future as DE becomes more 
dominant. It may provide an interaction between dark matter and DE and 
can drive transition from an early matter dominated era to a phase of ac- 
celerated expansion. Using the same phenomenon, DE and different forms 
of matter are discussed in the framework of F(T) theory which may help to 
discuss accelerated expansion of the universe. 

4.2 Using EoS Parameter 

Here we formulate some F(T) models in a slightly different way. We sub- 
stitute different values of parameter u in Eq. fl271) and solve it accordingly. 
Equation fl27|) can be written as 

AZF TT + [-4U(u + 1)+2E- AY] Ft + (to + 1)F = 0. (56) 

Now we construct F(T) models in the following cases. 
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Case 1: 

Putting lo — g in Eq.f l56|) . we have 

2ZF TT + (E — 2Y — ^U)F T + \f = 0. 

This has the following general solution 

" 8f/ _ 3 £ + 6 y + ^( 3jE ; _ gc7 - 6F) 2 - 48Z 

F(T) = c 3 exp T 



(57) 




+ C4 exp T 



where C3 and C4 are constants. 
Case 2: 

Consider the case when pressure is zero, i.e, u> = 0. In this case, Eq. (l56l) 
implies 

AZF TT + [-4U + 2E-4Y]F T + F = 0. (59) 
It has the following general solution 



F(T) = c 5 expT 



+ c 6 exp T 



6C/ - 3£ + QY + ^/{-6U + 3E - 6Y) 2 - 2AZ 



12Z 



QU -3E + QY - ^/{-6U + 3E - QY) 2 - 2AZ 



12Z 



(60) 



where C5 and Cq are arbitrary constants. 



Case 3: 

For oj — — 1, Eq.()5n]) becomes 

AZFtt + (2E - 4F)F T = 



(61) 



which has two solutions. The first solution yields F(T) = c-j while the second 
solution is given by 



F(T) = c 8 exp 



2Y -E 
2Z 



T 



(62) 
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where c 7 and c 8 are constants. Equations (155]) . ( 16"U|) and (152]) represent -F(T) 
models corresponding to radiation, matter and DE phases respectively. The 
exponential form of F(T) models represents a universe which always lies in 
phantom or non-phantom phase depending on parameters of the models [43J. 

5 Construction of EoS Parameter and Cos- 
mic Acceleration 

In this section, we derive EoS parameter by using two different F(T) models 
and also investigate cosmic acceleration. For this purpose, we evaluate p m 
and p m using the field equations and then construct the corresponding EoS 
parameter. 



5.1 The First Model 

Consider the following F(T) model 

F = aT+^, (63) 

where a and /3 are positive real constants. Replacing this value of F in 
Eqs.flU]) and JED, it follows that 

2k 2 Pm = (-AU + T)a + (3{l+AUT- l )T-\ (64) 
2K 2 p m = (2E -T)a- I3(1 + 2ET- 1 + 8ZT~ 2 )T~ 1 . (65) 

Dividing Eq. (l65|) by the EoS parameter is obtained as follows 

2(Y-U)(a-PT- 2 )-8ZpT-z 
(-4U + T)a + (3(1 + 4UT- 1 )T- 1 ' 1 ' 

Now we discuss this equation for particular values of a and j3. For a ^ 
0, — 0, we obtain 

w = + l ( i- W) ■ ( 67 ) 



3 V U . 

This leads to three different cases of u representing different phases of the 
evolution of the universe as follows: 
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• If jj > 1, then co < — 1 which corresponds to the phantom accelerating 
universe. 

• When jj < 1, the EoS parameter will be slightly greater than —1 which 
means that the universe stays in the quintessence region. 

• If jj = 1, we obtain a universe whose dynamics is dominated by cos- 
mological constant with u> = —1. 

It is interesting to mention here that model (1671) reduces to GR spatially flat 
Friedmann equation in the limiting case when anisotropy vanishes. The case 
a = 0, (3^0 does not provide meaningful results. 

5.2 The Second Model 

Assume F(T) has the form [22] 

F = aT + /3T n , (68) 

where n is a positive real number. The corresponding field equations will 
become 

2k 2 p m = (-4U + T)a + P(-AUnT~ l + l)T n , (69) 
2/t 2 pm = (2E - T)a + 2En(3T n ~ 1 - 4Zf3n(n - l)T n ~ 2 - f3T n . (70) 

Consequently, the EoS parameter takes the form 

2(U - Y)a + 2n/3(U - Y)T n ~ l + AZ(5n(n - l)T n " 2 

(_4c/ + T)a + ^(-^nT" 1 + l)T n ' ^ 

The case a ^ 0, /3 = leads to the same discussion as in the first case. For 
a = 0, (3 7^ 0, we have 



n 

w = -1 + 



2n + 1 



I <J-+ Z(;/ " l; 



(72) 



For any positive real number n, we can discuss as follows: 



When jj + Z ^ - > 1, Eq. (!72|) gives to < — 1 which represents the 
phantom accelerating universe. 
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For fj + 



Y , Z(n-1) 



U 2 



1, we obtain u = — 1 and hence the universe rests in 
DE era dominated by cosmological constant. 

• The case jj + z ^ 1 ^ < — 1 corresponds to the quintessence era because 
u > -1. 

Assuming n = 1 as a particular case in Eqs.(l69l) and ( JTOl) . we have 



Pm 
Pm 



— ( a + p)(-4U + T), 



2k 



2y a + (5)(2E-T). 



(73) 
(74) 



In the following, we discuss the evolution of the scale factor for Bianchi type 
I universe. For this purpose, we assume [22] 



Pm = + A (T) + Ax(r) Pr , 

Pm 



(75) 



such that A_i, Aq, A\ are constants. Substituting Eqs. (IT3l and (1741 in the 
above equation, it follows that 



2E-T 



-AU + T 



+ b + c(-AU + T), 



where 



(a + /3) 2 ' 
This equation leads to 

AU + 2E - 6 + 8C/c 



4k 4 A_x , 2k 2 /1 

6 = C = Ai. 



a + 



T 



± 



(76) 



(77) 



2(1 + c) 2(1 + c) 

x [b 2 - Aa - Aac - AEb + AE 2 + 8Ub - 1QEU + 16U 2 ] 1/2 . (78) 



Substituting this value of torsion in Eq. (j22j) . we have 

1 



H 



AU + 2E -b + 8Uc 
J 2(TTc) T 



2(1 + c) 



Vb 2 -Aa- Aac - AEb + AE 2 + 8Ub - WEU + 1QU 2 



1/2 



• (79) 
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The correspondingly average scale factor becomes 



R 



R Q exp { - 



AU + 2E - b + 8Uc 

J ; ; T 

2(1 + c) 



2 1 + c) 



x 



Vb 2 - 4a - 4ac - 4£6 + 4£ 2 + 8Ub - 1QEU + 1QU 2 



1/2 



dt k.80) 



As a special case of model ( 175|) . if we take A_i as a constant while A = 
= Al, we obtain standard Chaplygin gas EoS [44]. In this respect, Eqs.(|78 
and ( 1791) give the following results respectively 



T 
H 



(E + 2U) ± ^(E-2U) 2 
1 



J - (£ + 2C7) T ^{E-2U) 2 



The average scale factor for chaplygin gas has the form 
1 



R = R exp 



J — (E + 2U) =F vl^W 7 



(81) 
(82) 



(83) 



This represents an exponential expansion which may result a rapid increment 
between the distance of two non-accelerating observers as compared to the 
speed of light. As a result, both observers are unable to contact each other. 
Thus if our universe is forthcoming to a de Sitter universe [23J , then we would 
not be able to observe any galaxy other than our own Milky way system. 



6 Summary 

In this paper, we have investigated the recently developed F(T) gravity, 
where T is responsible for the cosmic acceleration without DE component. 
For this purpose, we have taken Bianchi type I spacetime which is one of the 
simplest models describing anisotropic, spatially homogenous and flat uni- 
verse. Some F(T) gravity models have been constructed by using two differ- 
ent approaches. In the first approach, we have used the continuity equation 
while in the second method, EoS parameter is used. These F(T) gravity 
models represent three different phases of the universe exhibiting different 
values of EoS parameter. The matter, radiation and DE eras respectively 
correspond to to — 0, u — 1/3 and to — — 1. 
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Matter dominated era describes expansion of the universe filled with non- 
interacting dust particles while radiation dominated era represents early uni- 
verse after the hot big bang. The DE case corresponds to the universe dom- 
inated by a strong negative pressure causing late-time acceleration. If we 
consider combination of radiation and matter, we may have more interesting 
results to study the developing universe. We have also constructed some 
models by using different combinations of EoS parameter. Also, we have 
obtained F(T) models in exponential form for some particular values of EoS 
parameter. 

Since the evolution of EoS parameter is one of the biggest efforts in obser- 
vational cosmology today. We have considered two well-known F(T) models 
and found the corresponding expression for u. These have been discussed for 
some particular values of parameters a, /3 which yield fruitful results corre- 
sponding to realistic situations. The cosmic acceleration has been discussed 
by using some F(T) models. We conclude that our universe would approach 
to de Sitter universe in the infinite future. It is interesting to mention here 
that one of the F(T) models (case 3) inherits a constant solution which may 
correspond to the cosmological constant. Notice that the isotropic expan- 
sion of the universe is obtained for A = which depends upon the values of 
unknown scale factors and parameters involved in the corresponding models 
[37]-|39j. It is worthwhile to point out here that our results correspond to 
FRW universe for the special case, i.e., A(t) = B(t) = C(t) = a(t). 
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